AFCRL  -  68-0249 


THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION 
OF  FOCAL  PLANE  ADAPTIVE  ANTENNA  TECHNIQUES 

BY 

MARVIN  COHN  AND  ROBERT  S.  LITTLEPAGE 


o 

CO 

00 


ADVANCED  TECHNOLOGY  CORPORATION 

1830  YORK  RD 

TIMONIUM,  MARYLAND  21093 


Contract  No.  F19628-67-C-0299 
Project  No.  4600,  8683 
Task  No.  460010,  868303 
Work  Uni».  No.  4600100’,  86830301 


FINAL  REPORT 

Period  Covered:  1  April  1967  to  31  March  1968 

April  1968 

Contract  Monitor:  Otho  E.  Kerr 
Microwave  Physics  Laboratory 

Distribution  of  this  document  is  unlimited.  It  may  be  released  to  the 
Clearinghouse  Department  of  Commerce,  for  sale  to  the  general  public. 


PREPARED 

FOR 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS  01730 


Reproduced  by  the 

CLEARINGHOUSE 

for  Federal  Scientific  &  Technical 
Information  Springfield  Va,  22151 


AFCRL  -  68-0249 


THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION  OF 
FOCAL  PLANE  ADAPTIVE  ANTENNA  TECHNIQUES 

by 

Marvin  Cohn  and  Robert  S.  Littlepage 


ADVANCED  TECHNOLOGY  CORPORATION 
1830  York  Road 
Timonium,  Maryland  21093 


Contract  No.  F19628-67-C-0299 
Project  No.  4600,  8683 
Task  No.  460010,  368303 
Work  Unit  No.  46001001,  86830301 


FINAL  REPORT 

Period  Covered:  1  April  1967  to  31  March  1968 

April  1968 

Contract  Monitor:  Otho  E.  Kerr 
Microwave  Physics  Laboratory 

Distribution  of  this  document  is  unlimited.  It  may  be  released  to  the 
Clearinghouse,  Department  of  Commerce,  for  sale  to  the  general  public. 


Prepared  for 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS 


ABSTRACT 


A  scalar  analysis  has  been  performed  to  determine  the  field 
structure  in  the  focal  region  of  apertures  focused  in  their  near 
field  (f/D  <  1).‘  Solutions  are  given  for  plane  waves'  at  normal 
incidence,  and  several  selected  angles  of  arrival.  The  results 
indicate  a  substantial  variation  from  the  Fraunhofer  far  field 
patterns.  Because  it  is  a  scalar  approach,  it  is  limited  to  large 
apertures  (D  >>\)  and  small  diffraction  angles. 

An  experiment  was  designed  to  complement  the  analysis. 

The  qualitative  agreement  is  very  good.  Sources  of  experimental 
error  have  been  found  which  account  for  the  quantitative  variations. 

A  technique  for  determining  the  field  structure  in  the  focal 
region,  due  to  a  reasonably  arbitrary  distribution  in  the  aperture 
plane,  is  given.  It  consists  of  the  superposition  of  an  ensemble  of 
plane  waves,  whose  amplitude,  relative  phase,  and  angle  of 
arrival  are  determined  by  the  aperture  plane  distribution. 
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1.  INTRODUCTION 

The  objective  of  this  effort  was  to  develop  a  manageable 
mathematical  expression  for  the  phase  and  amplitude  distribution 
in  the  focal  region  of  a  large  aperture  antenna  for  a  reasonably 
arbitrary  distribution  in  the  aperture  plane.  The  aperture  plane 
distribution  is  to  be  considered  the  result  of  either  reflector 
surface  errors  or  atmospheric  induced  errors.  Any  time  varying 
nature  of  the  aperture  plane  distribution  was  not  be  to  considered. 

The  result  of  the  increased  understanding  of  the  focal  region  fields 
can  be  used  to  improve  antenna  efficiency  by;  (  1)  either  the  design 
of  an  improved  single  feed  element  or  (  2)  an  array  of  feed 
elements. 

The  combined  effect  of  atmospheric  turbulence,  multipath, 
and  reflector  errors  results  in  extremely  complicated  equivalent 
aperture  distributions.  However,  regardless  of  the  complexity  of 
the  distribution,  a  Fourier  analysis  can  be  used  to  reduce  the 
problem  to  a  superposition  of  plane  waves  with  appropriate  phase 
and  amplitude.  Because  the  errors  are  in  general  not  periodic,  a 
Fourier  integral  would  be  required.  By  the  nature  of  atmospheric 
induced  errors,  and  for  reasonably  good  reflecting  surfaces,  most 
of  the  available  energy  will  be  confined  to  the  low  order  waves. 
Hence,  we  will  be  dealing  with  the  superposition  of  some  reasonable 
number  of  waves. 

Once  the  aperture  distribution  has  been  synthesized  by  an 
ensemble  of  plane  waves,  one  must  treat  the  diffraction  problem 
for  each  individual  wave.  This  amounts  to  finding  the  correct 
transform  function  for  the  particular  aperture. 

The  exact  solution  of  the  electromagnetic  boundary  value 
problem  for  diffraction  of  a  wave  at  an  aperture  in  an  opaque  screen, 


or  (  by  invoking  Babinot1  s  principle  of  complemontai  screens) 
diffraction  by  a  perfectly  conductirg  surface,  is  quite  tedious. 
Varying  degrees  of  approximation  have  been  applied  by  many 
investigators.  These  range  in  refinement  from  the  simple  scalar 
Fourier  transformation,  resulting  from  the  application  of  far  field 
approximations  in  the  Fresncl-Kirehhoff  integral,  to  the  vector 
treatment  of  Debye-Picht-Lunebcrg^  ^  with  integral  solutions  by 
the  method  of  stationary  phase. 

When  the  diffracting  aperture  has  focal  properties,  the 
problem  is  even  more  interesting.  This  is  the  case  that  is 
applicable  to  the  paraboloidal  and  spherical  reflecting  antennas. 
Until  recently,  most  antenna  designers  assumed  that  the  dif¬ 
fraction  pattern  in  the  focal  plane  of  the  dish  was  given  by  the 
Fourier  transformation  of  the  illumination  function.  In  other 

words,  the  field  structure  in  the  focal  plane  has  all  the  properties 

(2) 

of  the  far  field.  Sherman  has  shown  that  this  is  the  case  in  the 
Fresnel  region.  More  recent  analyses^  ^  and  experiments^  ^ 
have  indicated  that  this  is  not  the  case  for  apertures  focused  in 


1.  Kline,  M.  and  Kay,  I. ,  Electromagnetic  Theory  and  Geomet- 
rical  Optics,  Inter  science,  New  Y  u*k  (1963) ,  Chs.  XI  and  XII . 

2.  Sherman,  John  W. ,  "Properties  oi  Focused  Apertures  in  the 
Fresnel  Region,"  IRE  Trans,  on  Antennas  and  Propagation, 
Vol.  AP-10,  No.  4;  July  1962. 

3.  Hyde,  Geoffrey,  "Transverse  Focal  Region  Fields  of  a 
Spherical  Reflector,  "  Radio  Corporation  of  America, 

Scientific  Report  No.  2  on  AFCRL  Contract  No.  AF  19(628)- 
2758;  January  1966. 

4.  Loux,  P.C.,  Smilcn,  L.I.,  and  Martin,  R.W.,  "Transverse 
Antenna  Feed  Study,"  Sperry  Gyroscope  Company,  Final 
Report  on  AFCRL  Contract  No.  AF  19(  628)-2759;  January 
1966. 
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their  near  field  (small  focal  length  to  diameter  ratios) .  Since 
most  antennas  used  at  radio  frequencies  are  designed  with  --^■'s 
from  .25  to  .50,  it  is  particularly  important  that  the  diffracted 
fields  be  well  understood  in  order  to  design  proper  feeding 
structures. 

The  key  word  in  the  definition  of  this  study  is  "manageable". 
It  is  desired  to  develop  a  tool  for  designing  antenna  feed  systems. 
There  is  a  scale  associated  with  understanding  the  focal  plane 
fields.  Proceeding  from  the  Fraunhofer  approximation,  one  moves 
toward  an  exact  understanding  by  various  refinements  of  the 
analytical  technique.  In  terms  of  actual  feed  design,  this  scale 
has  a  point  of  diminishing  return.  Lacking  the  exact  analysis, 
the  merit  of  any  approach  rests  on  its  agreement  with  experiment. 

This  study  has  proceeded  on  the  basis  of  two  basic  restric¬ 
tions;  (l)  it  will  be  assumed  that  all  the  apertures  of  interest  are 
very  large  (D>>  \),  and  (2)  it  will  be  assumed  that  the  range  of 
diffraction  angles  of  interest  is  relatively  small.  It  is  difficult 
to  put  an  absolute  limit  on  the  latter  condition.  Since  the  goal  is 
to  develop  methods  of  designing  feed  structures,  the  extent  over 
which  physical  structures  can  be  located  will  be  limited  by  the 
incurred  aperture  blockage.  Hence,  knowledge  of  the  fields  far 
from  the  aperture  axis  is  not  required.  Another  advantage  of  the 
small  angle  limitation  is  that  most  of  the  available  energy  is  in  the 
principal  polarization.  A  vector  treatment  of  the  problem  may  be 
an  unnecessary  complication.  The  key  parameter  which  must  be 
retained  in  any  analytical  approach  is  the  effect  of  very  small  ^  . 
We  have  accomplished  this  by  a  more  rigorous  solution  of  the 
Fresnel-Kirchhoff  integral.  This  analysis  is  given  in  Section  2. 
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The  analysis  has  been  performed  for  circular  and  rectangular 
apertures.  Since  the  range  of  applicability  is  to  be  determined  by 
experiment,  it  was  felt  that  a  one  dimensional  aperture  offered  a 
simplification  for  the  initial  study.  For  this  reason,  an  exact 
solution  of  the  Fresnel-Kirchhoff  integral  equation  was  performed 
for  a  focused  infinite  slit.  This  geometry,  commonly  used  in  optics, 
has  an  excellent  long  wavelength  analog  in  the  parallel  plate -lens 
structure.  The  details  of  the  structure,  and  the  results  of  the 
measurements,  are  given  in  Section  3. 
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2.  ANALYSIS 


It  is  felt  that,  for  the  majority  of  antennas  of  interest,  much 
physical  insight  can  be  gained  by  a  more  rigorous  application  of 
the  scalar  diffraction  theory.  This  approach  neglects  polarization 
effects  and  hence  is  only  applicable  to;  (1)  the  principal  polarization 
in  the  region  of  the  geometric  axis  of  very  large  diameter  antennas, 
or  (2)  for  structures  which  by  design  permit  only  a  single  polari¬ 
zation.  The  analysis  is  applicable  to  a  very  large  class  of  antennas 
even  if  we  limit  the  aperture  size  to  fifty  wavelengths  or  greater 
(D/\  2:50).  It  has  been  shown^  ^  that  for  apertures  as  large  as 
this,  the  boundary  wave  due  to  edge  effects  has  a  negligible  contri¬ 
bution  at  moderate  diffraction  angles,  and  the  Fresnel-Kirchhoff 
integral  represents  an  excellent  approximation  to  the  diffracted 
fields.  For  a  two  dimensional  aperture  (three  dimensional  dif¬ 
fraction  problem),  this  integral  has  the  form, 


where 


us* 


*  ndS, 


(1) 


wave  disturbance  at  the  observation  point 
wave  disturbance  at  the  source  point 


ikr 

A 

—  —  s  Huygen  spherical  sources  in  the  aperture. 


5.  Wolf,  E.  and  Marchand,  E.W.,  "Comparison  of  the  Kirchhoff 
and  the  Rayleigh-Sommerfeld  Theories  of  Diffraction  at  an 
Aperture,"  J.  Opt.  Soc.  Am.,  54,  587  (1964). 

6.  Schivering,  F. ,  "On  the  Range  of  Validity  of  Fresnel-Kirchhoff  s 
Approximation  Formula,"  IRE  Trans,  on  Antennas  and  Propaga¬ 
tion,  Vol.  AP-10,  No.  1;  January  1962. 
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U  and  yU  *  n  are  taken  to  be  zero  everywhere  on  S  except  over 
b  b 

the  region  of  the  aperture.  Hence,  the  integral  over  the  closed 
surface  S  reduces  to  an  integral  over  the  aperture.  The  geometry 
associated  with  Equation  (1)  is  shown  in  Figure  1. 

Some  simplications  can  be  made  to  Equation  (1),  based  on 
the  size  of  the  aperture  but  irrespective  of  its  particular  geometry. 
Consider  the  term 


V 


ikr 
e _ 

r 


(2) 


For  the  apertures  of  interest,  one  can  define  a  range  of  possible 
values  for  the  diameter  (D),  focal  length  (f),  and  for  the  magnitude 
of  the  position  vector  r; 


50 

_f_ 

D 

r 


_D 

\ 

0.25 
f  . 


Considering  the  poorest  case,  one  still  finds  that  k  > >  Vr  by  an 
amount  6.283  >>  0.08.  Hence,  it  is  reasonable  to  neglect  the 
term  -Vr  in  Equation  (2). 

We  can  make  some  generalizations  about  the  form  of  Ug, 
The  focusing  properties  of  the  aperture  will  require  the  inclusion 
of  some  phase  distribution  (a  function  of  the  aperture  geometry) 
such  that  the  fields  due  to  every  Huygen  source  add  in  phase  at 


some  distance  f  along  the  aperture  axis,  for  a  wave  with  normal 
incidence  at  the  aperture.  The  incident  field  which  satisfies  these 
conditions  has  the  form 


E 

o 


ik*  r 
e 


+  ik8 


9 


(3) 


6 


where 


k  =  incident  wave  vector 
k5  =  phase  distribution  required  for  focusing. 
Equation  (1)  can  now  be  written  in  a  form  applicable  to  the  three 
dimensional  problems  of  interest; 


Eo  C  f  e*kr  p  ik*  r‘  +  ik6  ik  •  r'  f  ikS* 


ndS  . 


(4) 

It  should  be  noted  that  the  mechanism  for  including  the  focusing 
property  of  the  aperture  is  quite  arbitrary.  It  is  not  restricted  to 
any  particular  physical  structure  such  as  a  lens  or  shaped  re¬ 
flector.  Equation  (4)  has  been  studied  for  circular  and  rectangular 
apertures. 


2. 1  Solution  of  Equation  (4)  for  a  Circular  Aperture 

The  incident  wave  vector  k  is  shown  for  the  circular 
aperture  in  Figure  2 -a.  From  the  figure  it  is  apparent  that  wave 
vector  for  the  circular  aperture  has  components  of  the  form 


k  =  k  cos  iJj 
x 

k  =  -k  sin  i|*  cos  & 

y 

k  =  -k  sin  ^  sin  £ 
z 


(5) 

(6) 
(7) 


The  focusing  function  can  be  determined  by  considering  the  cross- 
sectional  view  of  the  aperture  shown  in  Figure  2-b.  Focusing 
requires  that  the  effective  propagation  paths,  and  hence  the  phase, 
from  every  point  in  the  aperture  to  the  point  z  =  f,  p'  =  o  be  the 
same.  This  requires  that 


where 


r  +  5  =  b  (constant) 


(8) 
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D  \2  2 

T)  tf.. 
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The  focusing  function  is  symmetric  in  the  variable  4>,  and  can  be 
found  from  Equation  (8); 


Equation  (11)  does  not  represent  any  particular  structure,  and  in 
fact,  when  one  applies  Equation  (11)  to  a  particular  case  care  must 
be  taken  in  defining  f  and  D, 

The  general  definition  of  the  geometric  parameters  of 
the  circular  aperture  are  given  in  Figure  3.  Under  the  coordinate 
transformations, 

*  =  *  (12) 

x  =  p  cos  $  (13) 

y  =  p  sin  $  ,  (14) 

the  form  of  the  incident  wave,  accounting  for  focusing,  becomes 

Tf  ~  i  k  [  z  cos  ^  -  P  ^  cos  ( <}>  -  £)  +  5  ( p)]  .  (15) 

US  "  o  6 

The  normal  component  of  del  operating  on  U  is  found  by  taking 


■us.  fl=[ 


where 


f  i 

p 

3Us 

-of  t  * 

3Us  cl  - 

17  kJ-“ 

(16) 

=  0, 

£•  n  =  -l 

• 

(17) 

If  the  aperture  plane  is  defined  to  be  the  plane  at  z  =  0,  then 


. 


i  k  [ -p  simjJ  cos  ( <J>  -  £)  H-  5  (p)] 


(18) 


V  U_  •  n  =  -  i  k  cos  i}j  E  c 
S  o 

Substitution  into  Equation  (4)  yields  the  field  due  to  a  circular 
aperture  of  diameter  D; 

D/2  2ir 

Up  =  Eq  j  f  [cos  ^  i-  cos  (n,  r)]  pdpdi  (19) 

p=0  4>=0 

where 

F(  p,<j>)  =  exp  [  i  k  (  r(p,<}>)  -  p  sinvji  cos  ($  -  C)  4  8  (p))  ]  (20) 


cos  (n,  r )  = 


r(p,<|>) 


(21) 


It  is  desired  to  determine  the  field  structure  in  the  focal  plane, 
defined  as  the  plane  perpendicular  to  the  z  axis  through  the 
point  z  =  f.  Hence,  from  Figure  3,  the  position  vector  r  has 
the  form, 


r  ( p. 4»)  =  [^f2  +  p2  +  p'2  -  2pp'  cos  (<(>  - , 


and  since 


one  can  write 


r(p,<»  = 


p'  =  f  tan  0 


f  2 

f  p  -  2  pf  tan  0  cos  (4>  -  <j>') 


2  . 
cos  0 


(22) 

(23) 


(24) 


Equation  (19) ,  with  Equations  (20)  ,  (21),  and  (24),  represents  the 
best  expression  for  the  focal  plane  field  structure  under  the  limita¬ 
tions  of  a  scalar  approach.  The  integral  cannot  bo  solved  in  closed 
form,  but  requires  a  numerical  computer  solution. 
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In  an  effort  to  simplify  Equation  (19)  in  order  to  gain 
physical  insight  into  the  problem,  some  further  approximations 
can  be  made.  Rewriting  Equation  (24)  in  the  form 

l 


r(p,4>)  =  —  +  pZ\ 

\Cos  0  • 


i  2  p  f  tan  9  cos  (§  -  ft1) 
f2  2 

- _  i  p 

cos  6 


.  (25) 


and  considering  only  small  0,  one  can  approximate  the  radical  by 
the  first  two  terms  in  the  binomial  expansion; 


r(p,ft) 


ter  * p2) 


p  f  tan  0  cos  (ft  -  ft1) 
f2  2 

- —  +  P 

cos  6 


(26) 


This  approximation  will  be  used  in  the  phase  term  (F(p.4>))  of 
Equation  (19).  The  amplitude  terms  (terms  other  than  F(p,ft)^  in 
Equation  (19)  can  be  approximated  by 


r(p,ft)  «  r(p)  =  (£2  t  p2) 


(27) 


since  they  are  slowly  varying  functions  of  0,  These  approximations 
result  in  an  integrand  which  is  a  separable  function  of  p  and  ft,  and 
Equation  (19)  takes  the  form 

D/2 


UP“-  4 


with 


ikEo  n  F(p)cik Ca(  pl  * 6< pD 


0  (f2tp2) 


COS  ft  + 


pdp  , 

(28) 


,  C  T.,/  .i  ,,  ,,  ,  pf  tan  0  cos  (ft  -  <f»')'\] 

F(p)  -  \  exp  -lk  fp  sin  ft  cos  (ft  -£)  *  ^ - y'J~~  1  dft 

0  L  /  f2  ,  .A  * 

- 2*  *  P  )  (i 

NCOS  0  / 
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=  Hr-2)5- 

\cos  0  / 


Equation  (  29)  can  be  solved,  and  the  solution  is  the  aero  order 


Bessel  function. 


where 


F(  p)  =  Ztt  Jq  n(  p)^ 


(p)  =  sin^  +  +  — - 


2  2  2  2  ~\  ■= 

P  i  tan  0  2p  f  tan  0  sin  ^  cos  ( t,  -  4> 1  )  | 

d  , .  ,  >r~  J- 


Hr  •  -■) 


The  resulting  expression  for  the  field  in  the  focal  plane  is, 


up*"T  Eo 


d/2  ikr(-4rt,2)*H 

U  L'  cos  0  J  p  -j 

\  “  J  (kn(p)]  cos  vpt - t  pdp  , 

o  (dtp2)  L  (dtP2fJ 


This  integration  still  requires  a  numerical  computer  solution.  How¬ 
ever,  depending  on  the  computer  capability,  the  integration  over  a 
single  variable  may  prove  to  be  a  valuable  simplification. 

2 . 2  Solution  of  Equation  (4)  for  a  Rectangular  Aperture 

The  geometry  for  the  rectangular  aperture  is  shown  in 
Figure  4-a.  The  incident  wave  vector,  k  ,  is  defined  in  Figure  4-b. 
By  an  anlysis  similar  to  that  for  the  circular  aperture,  one  can  show 
that  the  field  in  the  focal  plane  of  a  rectangular  aperture  with  sides 
a  and  b  is  given  by, 
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a/2  b/2 

U  3  .  lii  E  f  f  Ftx,y)  f 
P  4n  Eo  J  J  r,  2^FL 
-a/2  -b/2  If  +  x  +  y 


cos  4>'  cos  O'  t 


r—W] 

f2,  2  2  2J 

jf  fx  ty 


dxdy 


(34) 


where 


F(x,  y)  =  exp  i  k  ^r(x,  y)  -x  sin  <j>'  -y  cos  <}>'  sin  0*  +  8  (x,  y)^  . 

(35) 

For  this  geometry,  the  position  vector  and  focusing  function  are 
given  by. 


r(x,y)  = 


2  ,  2  . 
cos  9  cos  0 


-»(«£!*  run*) 


6(x,y)  = 


l 

2  2 

2 

a  +  b  .  ,2 

2,  2  2 

4 

x  t  y  +  f 

— 

J  (36) 
(37) 


Because  there  is  only  slight  interest  in  focused  rectangular  apertures, 
this  analysis  has  not  been  pursued  further. 

2.  3  Two  Dimensional  Diffraction  by  an  Infinite  Slit 

The  analysis  of  the  two  dimensional  problem  requires 
some  modification  of  the  Fresnel-Kirchhoff  equation  (1).  Assume 
a  volume  V  which  is  much  larger  in  length  than  in  any  radial 
dimension.  Hence,  cylindrical  coordinates  will  be  used  to  describe 
V.  The  pertinent  geometry  is  shown  in  Figure  5.  Assume  two 
functions,  U  and  v,  which  are  solutions  of  the  wave  equation  and 
hence  obey  Green's  second  identity; 


V2(U,v)  +  k2(U,v)  =  0 


(38) 


^  [v  VU  -  U  V  v  ]  •  ndS  =  0  . 


(39) 
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In  the  cylindrical  coordinates,  the  choice  for  v  which  satisfies  the 
wave  equation  is 


v  =  i  ir  H  (1)(kr)  , 
o 


(40) 


(1)  4 


where  H  is  the  zero  order  Hankcl  function  of  the  first  kind.  We 
o 

will  be  interested  in  the  asymptotic  forms, 


iirH  (1)(kr) 
o 


i(kr  -  17 /4) 


co 


(41) 


To  evaluate  the  disturbance  along  a  line  through  P,  of  length  L,  both 
asymptotic  regions  must  be  considered.  Therefore,  divide  the 
volume  V  with  a  sub-volume  defined  by  the  surface  S’.  Equation  (39) 
then  becomes 


§ 


[v  V  Ug-UgVv]  •  hdS 


^  [vVUs,-Us,Vv]  •  ndS* 


(42) 


Let  p  be  the  radius  of  the  cylinder  with  surface  S',  and  as  p  be¬ 
comes  small 


v 


i  it  H 


(1) 


(k  p)  -»  -2  In  (k  p)  . 


(43) 


The  normal  component  of  del  operating  on  v  results  in 


*  d  v 

vv- n  =  or 


(«) 


where  the  minus  sign  results  from  the  direction  of  n  being  out  of  the 
volume  V  as  shown  in  Figure  5.  Applying  the  Hankcl  identities  to 
Equation  (44) ,  along  with  the  small  argument  asymptotic  approxi¬ 
mation  yields 
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9  V  .  ,  T!  ^^1-  r  \  ->  - 

87=”klIl  (l'f,)  P  ■ 


Substituting  Equations  (43)  and  (45)  into  the  right  hand  side  of 
Equation  (42)  (which  is  the  integration  over  S'  where  the  above 
aoor oximations  arc  valid)  results  in 


vVU 


S'  S 


UqiV v]  *  ndS1  =  £  [2)n(kp)  VUg,]  •  ndS*  i  j  US‘(^)dS' 


In  the  limit  as  p 
only  one  term; 


becomes  vanishingly  small,  Equation  (46)  yields 


lim  <£>  U  ,(•-)  pd<ji  dz  ->•  4  tt  L  Ugt  , 
p-°  •1S,  p 

lim  <G  [ -2  ln(kp)vUg,]  •  n  pd«j»  dz  -  0  . 

p->°  Js, 

Therefore,  Equation  (42)  becomes, 

4  ir  LU  p  =  §  [vVUs-UsVv]  •  ndS  . 


Over  the 


3  surface  S,  the  large  argument  approximation  holds  for  v, 


(1)  ,  •  /itr"*  i(kr-u/4) 

=  i 11 H0  (kr)-*lYk?  G 


aw i 

-  =  .  Ti)  =  -iTTkH  (1)(kr)  cos(?,  n)  .  (51) 

Considering  the  large  argument  asymptotic  form  of  Equation  (51), 
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a  /2it  i(kr  -  3  ir/4)  .... 

Vvn  -*  -ik  i  /- —  e  cos  (v,  n)  . 

V  Kl* 


Substitute  Equations  ( 50)  and  (52)  into  Equation  (49)  to  show  that, 


up=  We 


i(kr-ir/4) 

•  e  [vU  .  n  -  ikUG  cos  (r,  h)]  r  d({>dz  . 

b  b 


Integration  with  respect  to  /.  causes  the  length  (L)  dependence  to 
drop  out  of  Equation  (53).  The  two  dimensional  equivalent  to  the 
three  dimensional  Frcsncl-Kirchhoff  integral  of  Equation  (1)  is, 

U  =  qi-  £  /pT  e  ^  1  ^{vUo.  n*ikU  cos  (r,  n)]  di  , 
P  4  ir  J  y  kr  » 


where 


rd<|>  ~  di  . 


The  line  i  is  the  perimeter  of  a  slice  of  the  volume  V  containing 

the  point  P.  As  in  the  three  dimensional  case,  we  will  assume  that 

U_  andv^c*  n  are  zero  everywhere  on  i  except  in  the  aperture. 

5  b 

This  reduces  the  integral  equation  (54)  to  an  integration  over  the 
aperture. 

A  cross-sectional  view  of  the  infinite  slit  is  given  in 
Figure  6.  The  incident  wave  will  be  of  the  form  given  in  Equation (3) 
and  this  particular  geometry  results  in 


U  =  E  e 
S  o 


ik  [  y  cos  O'  -  x  sin  O'  1-  6(x)  ] 


The  focusing  function,  5(x),  must  produce  a  spherical  phase  distri¬ 
bution  whose  radius  of  curvature  is  centered  at  the  point  y  ~  f,  x'  =  0 
Such  a  function  is  given  by, 
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A  more  convenient  form  of  Equation  (59)  can  be  generated  by 
normalizing  with  respect  to  wavelength; 


Equation  (64)  represents  the  complex  field  structure  of 
an  infinite  slit  focused  in  its  near  field.  The  integral  cannot  be 
solved  in  closed  form.  Since  it  must  be  applied  for  each  angle  of 
arrival  (O')  and  solved  for  a  high  density  of  values  of  diffraction 
angle  (0),  a  substantial  computer  effort  is  involved.  For  this 
reason,  a  perturbation  technique  has  been  investigated.  Expand 
Equation  (64)  in  a  Maclaurin  series  in  the  variable  0'  , 


V0')  =  Up(0'  =  0)  + 


yrfVl 

b  -  0G'n- 


O'  -  0 


The  first:  three  perturbation  terms  are, 
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This  approach  offers  the  £ amplification  that  integrals  are  solved  once 
for  some  range  of  diffraction  angle  0.  The  treatment  of  off  axis  waves 
is  then  reduced  to  the  algebraic  operation  of  Equation  (67).  This 
advantage  assumes  that  Equation  (67)  converges  for  some  reasonable 
number  k.  Since  the  various  ordered  derivatives  are  directly  multi¬ 
plied  by  0'  ,  this  condition  is  easier  to  satisfy  for  small  0'.  Some 
insight  into  the  rate  of  convergence  can  be  gained  by  considering 
Equation  (67)  at  the  particular  value  of  diffraction  angle,  0-0.  At 
this  value  of  0  the  integrands  of  the  various  perturbation  orders  have 
the  same  parity  as  the  order  number.  With  the  symmetric  limits  of 
integration,  wc  only  get  a  contribution  from  the  even  order  perturba¬ 
tion  terms.  This  yields  the  interesting  result  that,  in  order  to  have 
some  correction  at  all  values  of  0,  ki  2.  Another  result  of  using 
0  =  0  13  that  Ilia  integrands  are  nonoscillatory  functions,  and  hence 
solution  by  planimcter  measurement  is  quite  simple.  This  has 
been  done,  up  to  and  including  k  =  4,  for  an  f/D  =  .25  and  D/\  -  50. 

In  order  to  have  an  exact  value  for  comparison  to  the  perturbation 
approximation,  Equation  (64)  has  also  been  planimcter  integrated. 
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Al  (be  value  0  0,  the  only  contribution  to  this  integral  is  from  the 

cosine  term: 


(71) 

Although  this  is  an  oscillating  integrand,  and  requires  more  calcu¬ 
lated  points  to  define  the  curve,  for  small  O'  it  oscillates  very 
slowly.  Plots  of  the  integrand  for  some  small  values  of  O'  are 
given  in  Figure  7.  |lJp|2  from  Equation  (71)  has  been  compared 
to  the  perturbation  expansion  of  Equation  (67),  and  the  results  are 
given  in  the  following  table: 


llJpl2  (0-0) 


G'(°) 

Equation  (71  ) 

Equation  (67),  k  ~  2 

Equation  (67),  k  -  4 

0 

107.995 

— 

— 

.1 

105.808 

105.812 

105.820 

.5 

63. 339 

60.255  (.  2  dB) 

63.884  (.03  dB) 

1.0 

8.992 

.  01  36  (28.2  dB) 

12.717  (1.5  dB) 

1.5 

.211 

172.865  (29. 3  dB) 

29.113  (21.4  dB) 

The  numbers  in  parentheses  indicate  the  departure  of  the  perturba¬ 
tion  approximation  from  a  solution  of  the  exact  expression.  This 
shows  that  for  k  =  4,  the  perturbation  application  is  limited  to  angles 
less  than  1°.  To  apply  Equation  (67)  with  k  =  4  requires  a  computer 
solution  of  five  integral  equations.  In  light  of  the  limited  range  of 
O'  for  which  this  would  have  application,  it  was  decided  to  solve 
Equation  (64)  exactly  for  selected  values  of  O', 
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Equation  (64)  has  been  expanded  into  the  form, 
Up  s  A  I  i  B  . 

The  distribution  of  power  in  the  focal  plane  is  given  by, 
Power(0)(dB)=  10  Log1()|Up{0)|2  . 
and  the  phase  angle  is  given  by, 

Phase  Angle  (0)  =  tan  *  ‘X(q~)  * 


(72) 


(73) 


(74) 


A  computer  program  for  the  solution  of  Equations  (73)  and  (74)  is 
given  as  an  Appendix  to  this  report.  Computations  have  been  per¬ 
formed  for  diffraction  angles  ranging  from  -10°  to  +10°.  These 

were  performed  for  the  particular  cases  of  f/D  =  .25,  D/\  =  50, 

•  o  o  o  o  o  o 
and  0'=  0  ,  0.5  ,  1.0  ,  1.5  ,  2.0  ,  3.0  . 

The  solution  of  Equation  (73)  is  given  for  normal 

incidence  (O'  =  0)  in  Figure  8.  There  is  a  significant  difference 

between  this  pattern  and  the  usually  assumed  Fraunhofer  pattern. 

Within  the  range  of  0  for  which  computations  have  been  performed, 

there  is  a  spreading  of  the  pattern  and  substantially  higher  side- 

lobe  levels.  Figure  8  includes  a  scale  in  x/D  in  order  to  facilitate 

comparison  of  the  new  patterns  to  the  Fraunhofer  pattern. 

Interesting  points  of  comparison  are  summarized  in  the  following 


table: 

Fraunhofer  Approximation 

Equation  (73) 

3  dB  point 

0  =  0.56° 

0  =  0.90° 

1st  null 

0  =  1.15° 

0  =  2.0° 

1  st  sidelobc  level 

•13.48  dB 

-7.4  dB 

2nd  null 

0  a  2.3° 

0  =  4.6° 

2nd  sidelobc  level 

-17.92  dB 

-10.55  dB 

3rd  null 

0  =  3.45° 

0  =  7.25° 

3rd  sidclobe  level 

-20.84  dB 

-12.70  dB 

4th  null 

0  =  4.6° 

0  =  9.8° 
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FIG.8 -THEORETICAL  POWER  PATTERN  FOR  NORMAL  1NCID 


The  phase  distribution  in  the  focal  plane  is  given  for  this  case  in 
Figure  9.  It  is  very  similar  to  the  phase  behavior  of  the  Fraunhofer 
pattern,  with  discontinuous  steps  of  about  tt  radians  corresponding 
to  the  locations  of  nulls  of  the  power  pattern.  The  skirts  on  the 
leading  and  trailing  edges  of  the  phase  curves  are  probably  due  to 
computational  errors.  The  program  was  run  in  single  precision. 

In  the  vicinity  of  power  nulls,  the  functions  A  and  B  are  very  small 
numbers,  and  hence  computational  errors  are  more  significant. 

Very  interesting  comparisons  can  be  made  between  this 
solution  for  the  continuous  aperture  and  previous  results  for  an 
array  of  closely  spaced  elements.  The  continuous  infinite  slit 

(7) 

would  be  related  to  the  line  array.  Ricardi  has  analyzed  the 
case  of  a  line  array  focused  in  its  near  field,  with  emphasis  on 

determining  the  3  dB  spot  size.  Results  corresponding  to  our 

(81 

situation  (f/D  =  .25,  D/\  =  50)  are  given  by  Hansen'  '  and  indicate 
a  0.4\  spot  size  for  the  array.  This  would  correspond  to  a  dif¬ 
fraction  angle  of  .92°,  which  i6  in  excellent  agreement  with  our 
angle  of  .90°. 

The  power  patterns  for  O'  =  0.5°,  1.0°,  1,5°,  2.0°, 
and  3.0°  are  given  in  Figures  10  through  14.  Their  associated 
phase  distributions  are  given  in  Figures  15  through  19.  Two 
interesting  points  can  be  noted  from  the  patterns.  There  is  a 
gross  lack  of  symmetry  not  predicted  by  the  Fraunhofer  approxi¬ 
mation.  Also,  the  main  beam  is  shifted  by  an  angle  larger  than 

7.  Ricardi,  L.  J. ,  Electromagnetic  Theory  and  Antennas,  (E.C. 
Jordan,  ed.),  "Near-Field  Characteristics  of  a  Linear  Array", 
Proc,  URSI  Symp.  Electromagnetic  Theory  and  Antennas, 
Copenhagen,  June  1962,  Pergamon  Press. 

8.  Hansen,  R.C,,  Microwave  Scanning  Antennas,  Vol.  I, 
"Apertures",  Academic  Press,  New  York  (1964),  p,  44. 
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THEORETICAL  POWER  PATTERN  FOR  AN  ANGLE 
ARRIVAL  (2.0°) 


THEORETICAL  POWER  PATTERN  FOR  AN  ANGLE  OF 
ARRIVAL  (3.0°) 


PHASE  (RADIANS) 
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FIG.I7-THE0RETICAL  PHASE  FOR  ANGLE  OF  ARRIVAL  (1.5) 
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the  angle  of  arrival.  This  displacement  of  the  maximum  focal 
plane  field  intensity  is  1.6  times  the  angle  of  arrival.  The  phase 
distributions  still  contain  the  discontinuous  steps  of  about  tt 
radians,  with  a  superimposed  tilt  due  to  the  tilt  of  the  incident 
wave.  Further  comments  on  these  diffraction  patterns  will  be 
made  in  the  Experimental  and  Conclusions  sections  of  this  report. 
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3.  EXPERIMENT 

An  experiment  was  designed  to  measure  the  focal  region 
fields  of  an  aperture  focused  in  its  near  field.  A  structure 
equivalent  to  the  focused  infinite  slit  is  the  parallel  plate  guide 
containing  a  Ions.  This  structure  propagates  a  TEM  wave, 
polarized  perpendicular  to  the  parallel  plates,  and  only  permits 
diffraction  in  one  dimension.  These  polarization  and  diffraction 
properties  agree  with  the  analysis  of  Section  2. 

An  artist's  conception  of  the  structure  is  given  in  Figure  20. 

In  order  to  keep  the  physical  size  within  reasonable  limits,  the 
experiment  was  performed  at  35  GHz.  The  lens  is  piano- 
hyperbolic  and  is  constructed  of  Rexolite  1422,  having  an  index 
of  refraction  of  n  =  1.59  at  35  GHz.  The  spacing  of  the  parallel 
plates  is  .088",  which  insures  that  only  the  desired  polarization 
can  be  supported  even  in  the  region  of  the  lens.  The  total  height 
of  the  lens  is  16.87"  (50\)  and  the  height  of  the  parallel  plate 
structure  is  18.87".  The  two  inch  differential  in  height  is  evenly 
distributed  on  either  side  of  the  lens,  and  is  filled  with  absorber 
in  order  to  approximate  the  effect  of  the  parallel  plates  going  off 
to  infinity.  The  probe  consists  of  reduced  height  RG-96  waveguide 
to  fit  between  the  parallel  plates.  The  aperture  size  of  the  probe 
is  .030"  x  .080".  The  probe  is  equipped  with  a  drive  mechanism 
capable  of  accurate  movement  in  the  focal  plane. 

The  circuit  used  to  measure  the  diffraction  pattern  is  shown 
in  Figure  21.  The  source  is  a  1000  Hz  square  wave  modulated, 
klystron  oscillator.  The  phase  and  amplitude  of  the  diffracted  wave 
are. measured  by  a  bridge  circuit.  A  superheterodyne  receiver  is 
needed  in  order  to  operate  with  a  reasonably  small  probe  aperture. 
The  225'  spacing  between  the  source  and  the  lens  structure  more 
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DIELECTRIC  SPACER  COVERED 
WITH  ABSORBER 


FIG.20- PARALLEL  PLATE -LENS  ANTENNA 


FIG. 21-  SYSTEM  FOR  MEASURING  AMPLITUDE  AND  PHASE  OF  THE  DIFFRACTION  PATTERN 


than  satisfies  the  two  Rayleigh  distances  ( — - — )  required  for  plane 
wave  illumination.  Figure  22  is  a  view  of  the  pattern  range,  looking 
from  the  transmitter  toward  the  test  station.  The  actual  arrange¬ 
ment  of  the  structure  during  measurements  is  shown  in  Figure  23. 

In  Figure  24  the  parallel  plates  have  been  separated  in  order  to 
show  the  lens  and  probe. 

Initial  measurements  were  poor  due  to  several  effects.  Tests 
on  the  parallel  plate  structure  without  the  lens  indicated  the  presence 
of  a  large  standing  wave  due  to  reflections  from  the  rear  opening. 
This  was  corrected  by  the  addition  of  an  absorber  within  the  back¬ 
side  of  the  structure,  everywhere  except  where  probe  motion  was 
required.  Reflections  at  the  air-lens  interface  resulted  in  a  standing 
wave  within  the  lens.  This  produced  undesired  phase  variations  in 
the  diffracted  wave.  A  non-reflecting  coating  was  applied  to  the 
planar  surface  of  the  lens.  This  allowed  internally  reflected  energy 
to  pass  out  of  the  structure.  There  was  also  a  problem  with  the 
lens  properly  contacting  the  parallel  conducting  plates.  This  was 
overcome  by  chemically  depositing  a  thin  film  of  silver  on  the 
broad-flat  surfaces  of  the  lens. 

Patterns  were  measured  for  normal  incidence  and  for  all  the 
angles  of  arrival  for  which  computations  were  performed.  The 
results  are  given  in  Figures  25  through  30.  (In  some  cases,  only 
the  maximum  of  a  sidelobe  could  be  measured,  and  this  is  shown 
as  a  dot.)  A  point  by  point  comparison  with  the  calculated  data 
indicates  rather  poor  agreement,  although  always  better  than  that 
obtained  by  comparison  with  calculations  based  on  the  Fraunhofer 
approximation.  However,  the  qualitative  agreement  is  very  good. 
This  is  particularly  evident  for  the  highly  unsymmetric  pattern  for 
the  angle  of  arrival  O'  =  3.0°.  Two  other  interesting  comparisons 
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FIG.  22  -  A  VIE/V  OF  THE  PATTERN  RANGE 
LOOKING  FROM  THE  TRANSMITTER  SITE 
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FIG.  23  -  A  VIEW  OF  THE  MEASURING  APPARATUS 


FIG.  24  -  THE  TEST  STRUCTURE  WITH  THE  PLATES  SEPARATED  TO  SHOW  THE  LENS  AND  PROBE 
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MEASURED  POV/ER  PATTERN  FOR  AN  ANGLE  OF  ARRIVAL 


EASURED  POWER  PATTERN  FOR  AN  ANGLE  OF  ARRIVAL 


N  FOR  AN  ANGLE  OF  ARRIVAL 


MEASURED  POWER  PATTERN  FOR  AN  ANGLE  OF  ARRIVA 


The 


can  be  made  between  the  calculated  arid  measured  patterns, 
calculated  patterns  were  presented  in  a  normalized  form.  The 
computer  solution  was  in  a  non-normalized  form  permitting  direct 
comparison  of  power  for  different  angles  of  arrival.  A  curve  of 
the  power  at  the  maximum  of  the  diffraction  pattern  (relative  to 
that  at  0  "  0  when  O'  =  0)  versus  the  angle  of  arrival  is  given  in 
Figure  31.  The  measured  data  is  in  very  good  agreement  in  both 
magnitude  and  shape.  The  other  comparison  is  the  displacement 
of  the  angle  of  the  diffraction  maximum  (  0  MAX)  versus  the  angle 
of  arrival.  The  Fraunhofer  approximation  would  predict  a  one-to- 
one  relationship.  Our  analysis  predicts  that  for  angles  of  arrival 
up  to  about  2.0°,  the  ratio  of  displacement  is  1.6.  The  measured 
data  yields  a  ratio  of  2.2.  These  results  are  shown  in  Figure  32. 
Measured  phase  data  has  not  been  included  in  these  figures  because 
the  system  eignal-to-noisc  ratio  did  not  permit  accurate  phase 
measurement.  However,  the  abrupt  steps  of  about  ir  radians  were 
observed  to  correspond  to  the  position  of  power  nulls. 

There  is  more  than  a  casual  agreement  between  this  measured 
data  and  the  calculations.  There  arc  two  fundamental  departures 
from  the  calculations;  (1)  the  measured  patterns  are  spread  in 
comparison  to  the  calculations,  and  (2)  the  angular  displacement 
of  the  diffraction  maxima  are  too  large.  This  hus  prompted  an  in 
depth  study  of  the  experiment  for  sources  of  error,  their  potential 
magnitude,  and  the  direction  of  the  resulting  effects. 

As  with  any  experiment  involving  field  probing,  the  probe 
must  be  considered  as  a  source  of  experimental  error.  Because 
of  our  limited  receiver  sensitivity,  we  were  forced  to  operate  with 
a  probe  aperture  which  subtended  an  angle  of  1.1°  in  the  focal  plane. 
The  integrating  effect  of  the  probe  makes  the  accurate  measure¬ 
ment  of  the  position  of  nulls  and  the  particular  shape  of  the  pattern 
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ANGULAR  DISPLACEMENT  OF  THE  MAIN  BEAM  VS 
ANGLE  OF  ARRIVAL 


(  3  dB  points,  etc.)  very  difficult.  Also,  the  manual  movement  of 
the  probe  necessitated  a  quantized  measurement  of  the  pattern. 

The  data  was  presented  exactly  as  measured,  and  hence,  some  of 
the  sidelobes  appear  pointed. 

An  even  more  fundamental  problem  exists  in  the  design  of 
the  lens.  In  an  effort  to  avoid  amplitude  tapering  due  to  loss  in 
the  lens  material  (the  analysis  assumed  uniform  amplitude 
illumination) ,  a  low  dielectric  constant  material  was  used.  This 
resulted  in  a  thick  lens,  for  which  the  principal  plane  (effective 
aperture  plane)  location  is  not  accurately  known.  The  indications 
are,  from  applying  techniques  of  ray  optics,  that  the  effective 
aperture  of  our  lens  is  less  than  the  desired  50X. 

A  further  complication  arises  from  the  radiation  pattern  of 
the  probe  itself.  An  isolated  measurement  of  the  probe  revealed 
that  it  introduces  a  3  dB  taper  out  to  the  45°  point  of  its  radiation 
pattern.  This  superimposes  an  amplitude  taper  across  the  aperture 
and  would  have  the  effect  of  spreading  the  pattern. 

Because  of  the  mutual  interplay  between  the  probe  and  the 
lens,  it  was  felt  that  the  question  of  determining  the  effective 
aperture  could  best  be  answered  by  direct  measurement  of  the 
entire  structure1  s  far  field  pattern.  With  the  probe  fixed  on  the 
axis  of  the  lens,  rotation  of  the  entire  structure  should  produce  the 
Fraunhofer  far  field  pattern.  The  results  are  given  in  Figure  33. 
The  first  null  occurring  at  about  2°  indicates  that  the  effective 
aperture  is  about  30  X.  The  location  of  the  focal  plane  of  the 
structure  was  verified  by  measuring  a  series  of  diffraction  pat* 
terns  through  several  planes,  and  determining  the  plane  for  which 
the  width  of  the  main  beam  was  a  minimum.  This  corresponded 
exactly  to  the  predicted  location  of  the  focal  plane.  The  combined 
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result  is  that  the  experimental  structure  has  an  i/D  too  large  and 
a  D/k  too  small.  The  i/D  is  approximately  .  57  which  is  still  in  the 
near  field  and  hence  in  the  region  of  interest.  The  D/\  is  approxi¬ 
mately  30.  Lacking  a  solution  of  Equation  (64)  for  this  value,  we 
can  only  say  that  it  is  in  the  correct  direction  to  account  for  the 
variations  in  the  measured  and  calculated  patterns. 

This  problem  was  alluded  to  in  Section  2  when  it  was  noted 
that  caution  is  required  in  going  from  the  mathematical  model  of 
a  focused  aperture  to  some  real  structure.  The  model  assumed 
that  all  phasing  was  accomplished  in  an  infinetesimal  region 
called  the  aperture  plane.  This  is  not  the  case  in  lenses  or  re¬ 
flecting  antennas.  The  application  of  Equation  (64)  to  some 
particular  structure  will  first  involve  determining  the  correct  f 
and  D.  Another  result  of  the  focusing  model  is  its  inability  to 
correctly  account  for  the  skewed  path  through  a  lens  or  parabola 
due  to  off  axis  Ulumination.  To  include  this  effect  would  amount 
,o  an  unnecessary  complication  for  the  small  angles  of  arrival 

which  are  of  interest. 


4.  CONCLUSIONS 


Let  us  assume  in  this  discussion  that  all  the  apertures  of 
interest  are  very  large  (D>>\.)  and  that  the  diffraction  angles  of 
interest  are  small.  As  previously  stated,  this  allows  one  to  neglect 
polarization  (essentially  treating  the  principal  polarization)  and  con¬ 
fine  ail  further  discussion  to  the  scalar  approximation.  Under  this 
assumption,  the  exact  solution  for  the  focal  region  fields,  resulting 
from  complex  amplitude  and  phase  distributions  in  the  aperture 
plane,  requires  solving  an  integral  of  the  product  of  the  aperture 
plane  distribution  and  the  aperture  transformation  function.  This 
requires  a  computer  solution  for  each  aperture  plane  distribution. 

We  feel  that  a  more  manageable  approach  is  the  statistical  combina¬ 
tion  of  an  ensemble  of  plane  waves.  The  weighting  factor  (the 
relative  amplitude  and  phase  between  component  waves)  is  determined 
by  the  nature  of  the  aperture  distribution.  In  principal,  the  com¬ 
ponent  plane  waves  could  have  any  angle  of  arrival  from  the  forward 
hemisphere.  This  will  also  be  a  function  of  the  type  of  aperture 
error.  The  simplification  of  this  approach  arises  from  the  fact  that 
once  the  integral  of  the  aperture  transform  function  is  solved  for  a 
given  angle  of  arrival,  different  aperture  plane  distributions  can  be 
solved  by  an  algebraic  manipulation  of  the  component  waves, 
weighted  with  appropriate  phase  and  amplitude. 

Prior  to  treating  particular  aperture  plane  distributions,  it 
was  necessary  to  develop  the  aperture  transform  function  for 
apertures  focused  in  their  near  field.  Very  significant  results  came 
from  this  investigation.  The  diffraction  characteristics  of  a  single 
plane  wave  were  substantially  different  from  the  Fraunhofer  approxi¬ 
mation.  The  main  beam  is  much  wider  and  the  sidelobe  level  is 
higher  than  customarily  assumed.  For  very  small  angles  of  arrival, 


the  shape  of  the  main  beam  remains  constant,  but  it  is  displaced 
by  an  angle  greater  than  the  angle  of  arrival.  This  is  significant 
to  the  design  of  an  adaptive  feed  system  since,  even  through  the 
aperture  plane  errors  only  excite  small  angular  waves,  these 
angles  are  exaggerated  by  diffraction.  This  could  result  in  the 
angle  subtended  by  the  feed  structure  being  considerably  larger 
than  the  cone  of  angles  of  the  plane  waves  required  to  synthesize 
the  aperture  distribution.  The  final  design  of  the  feed  system  will 
be  predicated  on  the  granularity  of  the  interference  pattern  of  this 
ensemble  of  plane  waves.  This  will  determine  the  number  of 
elements  required  to  recover  a  given  percentage  of  the  energy  in 
the  focal  plane.  The  size  and  location  of  the  elements  in  the  focal 
plane  must  be  chosen  so  as  to  produce  the  minimum  disturbance 
of  the  focal  region  fields. 

The  final  determination  of  the  focal  region  field  structure 
for  a  particular  aperture  distribution  will  require  the  complex 
addition  of  the  constituent  plane  waves.  There  has  not  been 
sufficient  time  to  investigate  these  ensemble  characteristics 
during  this  contract.  However,  one  particular  aperture  distri¬ 
bution  of  interest  has  been  set  up  to  demonstrate  the  technique. 

A  periodic  aperture  distribution,  with  a  period  greater  than  the 
v/avelength  of  the  incident  radiation,  can  be  represented  by  the 
superposition  of  two  plane  waves.  The  magnitude  of  the  aperture 
plane  variations  will  be  determined  by  the  relative  amplitude  of 
the  two  incident  plane  waves.  The  period  of  the  aperture  plane 
variations  will  be  determined  by  the  angles  of  arrival.  Assume  a 
reference  wave,  which  after  diffraction  has  the  form 

UP  <V)  s  Ao(0i,)  +  iBO(ei,)  •  (75) 

o 
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The  second  wave,  which  produces  the  desired  interference  pattern, 
has  the  form, 


<V)tiBl(e2 


i  A<|> 


(76) 


where  is  the  initial  phase  difference  between  the  two  waves. 

This  determines  where  in  the  aperture  plane  the  maxima  of  the 
periodic  variations  occur.  The  magnitude  of  the  variations  will 
depend  on  the  amplitude  ratio  a  ,  where, 

t  1 

Ei<Y>  *  ‘iW  :  <77> 

resulting  in, 

W1"*  "lVY1  ■  W^lW  ’  (78) 

The  A' s  and  B's  are  the  real  and  imaginary  solutions  of  Equation(64) 
for  the  particular  values  of  0'.  The  total  resulting  field  in  the  focal 
plane  is  given  by. 


Up  =  A  4-  i  B  , 


with 


A  =  A 


B  =  B 


^Oj')  +  a^A^O^')  cos  "  Bo^02^  8in 

,(0^)  +  aL  [^Ao(02')  sin  A<|>  +  BJ0.,')  cos 


(79) 

(80) 
(81) 


The  diffraction  pattern  can  be  found  from  Equations  (80)  and  (81), 
and  is  given  by, 

|opr  =  K>  (9,')|2  +  -^lup  (92')|2  *  2a,  (A0(e,')fAo(e2')  co,  A  + 
o  o  V.  *- 

-  BQ(®2^  8*n  *  Bo^i'm  a0^2 s*n  *  ^0^2'^  cos  }• 

J  L  (82) 
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The  zero  subscripts  refer  to  the  calculations  performed  in  Section  2 
assuming  unit  incident  amplitude,  E^  =  1,  The  first  two  terms  on 
the  right  hand  side  of  Equation  (82)  are  ordinary  diffraction  patterns 
for  particular  angles  of  arrival  (0*)»  like  those  presented  in  Fig¬ 
ures  10  through  14.  The  last  term  in  Equation  (82)  would  require 
some  algebraic  calculations.  The  necessary  A's  and  B's  have  been 
included  as  print  outs  in  the  computer  program. 

Even  this  simple  example  demonstrates  that  for  combining 
more  than  two  waves,  the  algebra  can  become  quite  tedious.  How¬ 
ever,  in  principal  it  can  be  solved  on  a  desk  calculator,  whereas 
the  exact  problem  cannot.  In  some  cases,  when  many  waves  are  to 
be  investigated,  it  would  be  simpler  to  include  the  performance  of 
Equation  (82)  in  the  computer  program. 

Aside  from  developing  the  tool  for  treating  the  problem  of 
distorted  apertures,  the  analysis  in  this  report  provides  an 
excellent  approximation  for  the  behavior  of  assumed  perfect 
apertures  focused  in  their  near  field.  The  analysis  was  tested  by 
performing  experiments,  and,  although  the  experimental  structure 
did  not  exactly  correspond  to  the  parameters  used  in  the  analysis, 
the  qualitative  agreement  was  good.  The  analysis  predicted  that 
for  an  f/D  =  .25  and  D/\  =  50,  the  main  beam  would  be  1.74  times 
larger  than  the  Fraunhofer  approximation.  Other  experimenters 
have  noted  and  reported  this  effect.  Loux  et.  al.  measured 
patterns  in  the  near  field  of  a  parabolic  dish  with  f/D«  •  33  and 
D/\  a  50.  They  found  that  the  radius  of  the  Airy  disc  was  1.5 
times  the  predicted  (Fraunhofer)  value.  As  f/D  increases,  our 


analysis  yields  results  which  approach  those  predicted  by  the 
Fraunhofer  approximation.  Hence,  the  variation  in  the  1.5  ratio 
of  Loux  and  our  1.74  is  accounted  for  by  the  f/D  going  from  .33  to 
.25.  The  enlargement  of  the  main  beam  and  the  increased  sidelobe 
level  are  sufficient  to  justify  attempting  to  recover  this  energy. 
Based  on  the  over  simplified  approach  of  planimetering  the  area 
under  a  focal  plane  curve  of  absolute  power,  it  appears  that  a  feed 
designed  to  properly  add  the  main  lobe  and  the  first  sidelobes 
would  have  an  efficiency  gain  of  2  -  2.  5  dB  over  a  single  feed  de¬ 
signed  to  collect  the  main  lobe  predicted  by  the  Fraunhofer  approxi¬ 
mation.  Because  the  off  axis  illumination  indicated  highly  un- 
symmetric  sidelobe  levels  and  a  steering  factor  of  1.6,  these 
considerations  should  yield  improvements  in  the  design  of  monopulse 
and  conical  scan  systems. 

A  potentially  interesting  application  of  the  results  of  the 
analysis  is  its  use  in  the  design  of  feed  structures  to  synthesize 
prescribed  far  field  radiation  patterns.  Given  a  pattern  require¬ 
ment,  one  can  approximate  it  as  an  ensemble  of  incident  waves 
whose  individual  amplitudes  and  phases  correspond  to  those  of  the 
required  pattern  at  the  appropriate  angles  of  arrival.  Using  the 
aperture  transform  developed  during  this  program,  the  focal  plane 
field  distributions  of  each  of  these  incident  waves  are  calculated 
and  superimposed.  A  feed  assembly  is  then  designed  to  approxi¬ 
mate  this  focal  plane  field  distribution.  This  is  analogous  to  the 

(10) 

work  of  Woodward  for  pattern  synthesis  by  arrays. 


10.  Woodward,  P.M. ,  "A  Method  of  Calculating  the  Field  Over  A 
Plane  Aperture  to  Produce  a  Given  Polar  Diagram,"  Journal 
of  the  Institute  of  Electrical  Engineers,  93,  Part  III  A,  1947, 
pp.  1554-1558. 
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5.  RECOMMENDATIONS 

The  aperture  transform  function  developed  during  this  study 
provides  an  excellent  approximation  to  the  focal  region  fields  of 
apertures  focused  in  their  near  field.  It  provides  the  basic  tool  for 
the  study  of  focal  region  field  dependence  on  complex  aperture 
plane  distributions.  Although  time  did  not  permit  an  in  depth  study 
of  aperture  distributions,  the  effort  has  progressed  to  the  point 
where  data  pertinent  to  focal  plane  adaptive  combining  can  now  be 
generated.  A  program,  which  makes  maximum  use  of  what  has 
been  developed,  will  be  given  for  the  continuation  of  study. 

Since  the  antennas  of  interest  are  of  the  large  dish  variety, 
computations  should  be  performed  for  Equation  (33)  of  the  circular 
aperture.  This  will  determine  the  nature  of  sidelobes,  steering 
angle  ratio,  and  phase  relationships,  similar  to  the  data  presented 
for  the  focused  infinite -slit.  This  type  of  data  is  necessary  in 
order  to  determine  the  resulting  field  granularity  from  the  super¬ 
position  of  the  several  waves  required  to  synthesize  the  aperture 
plane  distribution.  The  design  of  an  array  of  feed  elements  is 
based  on  the  maximum  recovery  of  energy,  while  introducing  a 
minimum  disturbance  to  the  fscal  region  fields.  This  may  not 
correspond  to  a  plane  array  located  in  the  focal  plane  of  the  dish. 
The  optimum  size  and  location  of  elements  can  only  be  determined 
when  the  fields  in  the  focal  "region”  are  known.  It  is  felt  that  a 
three  dimensional  picture  of  the  field  structure  is  needed.  This 
can  he  accomplished  by  the  simple  adjustment  of  Equation  (33)  to 
solve  for  the  field  structure  in  cuts  other  than  the  focal  plane. 

This  would  also  answer  the  question  of  how  critical  is  the  location 

(2) 

of  individual  elements.  It  is  known  that  the  "depth  of  focus" 
(distance  between  the  axial  3  dB  points)  decreases  with  decreasing 
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focal  length  in  the  Fraunhofer  approximation.  This  situation  may 
be  more  acute  for  apertures  focused  in  their  near  field  since  both 
amplitude  and  phase  determine  the  diffraction  pattern.  Some 
further  advantage,  in  terms  of  aperture  blockage,  might  be 
realized  by  locating  the  elements  over  a  curved  surface.”  Again, 
Equation  (33)  is  easily  set  up  to  solve  for  this  case,  provided  the 
curved  surface  is  reasonably  regular. 

Since  it  is  known  that  the  Fraunhofer  approximation  for  the 
diffracted  fields  is  good  for  apertures  focused  in  the  Fresnel  region, 
it  would  be  desirable  to  have  a  better  determination  of  where  the 
Fresnel-near  field  boundary  occurs.  This  can  be  accomplished  by 
solving  for  other  f/D  ratios.  This  would  also  make  the  results  of 
this  analysis  applicable  to  a  wider  variety  of  dishes. 

In  many  cases,  the  antenna  designer  is  not  concerned  with 
aperture  plane  errors.  This  would  be  the  case  where  the  reflector 
surface  is  reasonably  good  at  the  operating  frequency  and  the 
aperture  size  is  small  enough  that  atmospheric  induced  errors  are 
not  substantial.  It  might  be  desired  to  amplitude  taper  in  order  to 
reduce  sidelobe  level  or  alter  the  main  beam  characteristics.  Our 
analysis  provides  a  more  accurate  aperture  transform  function 
than  the  Fourier  transformation.  The  unit  amplitude  distribution, 
assumed  in  our  analysis,  can  be  replaced  by  certain  simple  distri¬ 
butions  (such  as  the  cosine  on  a  pedestal)  without  appreciably 
changing  the  computer  effort  to  solve  Equation  (33).  This  would 
determine  the  phase  and  amplitude  distribution  required  in  the 
focal  plane  to  produce  the  desired  aperture  plane  amplitude  distri¬ 
bution. 

It  will  do  desired  to  experimentally  verify  many  of  the  above 
calculations.  During  the  experiments  already  performed,  it  was 


found  that  the  existing  thick  lens  was  a  poor  approximation  to  the 
model  for  which  calculations  were  performed.  Since  the  structure 
represents  a  simple  technique  for  testing  the  limits  of  applicability 
of  (he  analysis,  it  is  recommended  that  the  existing  lens  be  re¬ 
placed  by  a  thin  lens,  whose  effective  aperture  and  focal  length  are 
more  easily  determined.  The  geometry  will  remain  plano- 
hyperbollc  because  of  the  reduced  machining  when  only  one  curved 
surface  is  involved.  The  use  of  a  higher  permittivity  material  will 
reduce  the  lens  thickness  to  a  more  reasonable  size.  As  an  example, 
Emerson  &  Curning,  Inc.,  Hi  K  #140  has  a  dielectric  constant  of  20 
and  a  loss  tangent  of  .0009.  A  lens  with  f/D  =  .25  and  D/\  =  50 
would  have  a  maximum  thickness  of  3.89X.  The  amplitude  taper 
introduced  by  lens  loss  would  be  0.43  dB. 

Based  on  the  results  of  this  program,  the  use  of  a  corrected 
lens  should  result  in  very  good  agreement  between  the  analysis  and 
experiment.  Once  the  fundamentals  have  been  verified  for  single 
waves  at  various  angles  of  arrival,  one  might  tend  to  assume  that 
the  technique  of  superposition  is  so  well  known,  that  further  experi¬ 
mental  verification  is  not  required.  This  may  not  be  the  case.  The 
scalar  analysis  neglects  certain  physical  phenomena  on  the  basis  of 
negligible  contribution  for  large  apertures  and  small  diffraction 
angles.  Noteworthy  is  the  surface  wave  generated  in  the  plane  of 
the  aperture.  We  know  that  the  excitation  efficiency  of  this  wave 
increases  with  increasing  angle  of  arrival.  Effects  which  were 
insignificant  for  a  single  wave  may  be  enhanced  by  the  complex 
addition  with  other  incident  waves.  For  this  reason,  it  is  felt  that 
the  focal  region  fields  should  be  measured  for  some  small  number 
(2  or  3)  of  incident  waves. 
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One  problem,  which  is  very  difficult  to  treat  analytically, 
but  which  could  be  exhaustively  studied  experimentally,  is  the 
effect  of  aperture  blocking  by  the  feed  assembly.  The  use  of  a 
transmission  structure,  as  the  parallel  plate  lens,  has  allowed 
the  probing  of  the  diffracted  fields  without  any  blockage  effect. 
However,  controlled  amounts  of  blockage  could  be  introduced  in 
front  of  the  lens,  and  the  effect  determined  by  field  probing.  The 
trends  in  the  perturbation  of  the  diffracted  fields  would  be  measured 
as  a  function  of  the  size  and  location  of  the  blockage.  Since  the 
effect  of  the  blockage  is  to  produce  errors  in  the  aperture  plane 
distribution,  these  measurements  would  also  provide  an  indication 
of  how  fast  the  focal  region  fields  degrade  with  aperture  error.  It 
would  provide  insight  as  to  the  number  of  plane  waves  required  to 
reasonably  approximate  typical  aperture  plane  errors. 

It  is  felt  that  there  is  now  available  sufficient  data  to  design 
an  improved  feed  for  a  single  on-axis  wave.  This  could  potentially 
increase  the  antenna  efficiency  by  a  couple  of  dB.  A  small  array 
of  properly  designed  near  field  feeds  could  be  used  to  reduce  the 
sidelobe  level.  The  resulting  reduction  in  antenna  noise  for  low 
noise  systems  could  result  in  benefits  *hat  '•xceed  those  predicted 
by  the  increase  in  antenna  efficiency. 
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APPENDIX 


C  DIFFRACTION  BY  A  LINEAR  APERTURE  F/D=.  25, 
1D/LAMBDA=50 

3001  FORMAT(  15X,F4.1,9X,2{  F25. 8) ,  F23.4,  F25.4) 

3002  ^ORMAT(*l*,  15X,  *THETA*,  24X 
1  POWER*,  20X,  *  PHASE*) 

PRINT  3002 


THETA=-10. 1 
DX=.  005 
F=12. 5 
D=25. 

THETAP=3.  0 
TP=THETAP*.  0174523925 
SP=SINF(TP) 

CP=COSF(TP) 

1=0 

PI=3. 1415926536 
PI4=Pl/4. 

FSQ=F**2 
PI2=2.*PI 

DEL=SQRTF(D**2+FSQ) 

20  1=141 
A=0. 

B=0. 

THETA=THETA+.  1 
T=THETA*.  0174523925 
SN=SINF(T) 

CS=COSF(T) 

TAN=SN/CS 
G1=FSQ/(CS**2) 

FTAN=F*TAN 
X=-D 
XSQ=X**2 
DELTA=DEL-SQRTF(XSQ4FSQ)-X*SP 
G=SQRTF(GltXSQ-2.*X*FTAN) 
Z=CP+F/G 

TARG=P12*( G+DELTA)  -PI4 
S=SINF(TARG) 

C=COSF(TARG) 

GRT=SQRTF(G) 

AGRND=C*Z/GRT 
BGRND=S*Z/GRT 
DO  100  J  =  1,5000 
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A*,  24X, *B*,  24X, * 


.  ! 

i 

*  - 


. 


J'  - '  . 

- - - — ^ 


X=X+.01 

XSQ=X**2 

DELTA-DEL-SQRTF(XSQ+FSQ)-X*SP 
G=SQRTF(GHXSQ-2.  *X*FTAN) 

Z=CPf f/g 

'  TARG=PI2«(G+DELTA)-PI4 
S=SINF(TARG) 

C=COSF(TARG) 

AGRNDO=AGRND 

BGRNDO=BGRND 

GRT=SQRTF(G) 

AGRND=C*Z/GRT 

bgrnd=s*z/grt 

A=A+DX*(AGRND+AGRNDO) 

100  B=B+DX*(BGRNDf  BGRNDO) 

USQ=.  25*(A#*2+B**2) 

POWER-10.  ^ALiOGlO(USQ) 

PHASE=ATANF(B/A) 

PRINT  3001,  THETA,  A,  B,  POWER,  PHASE 
IF  (1-200)20,  20, 140 
140  RETURN 
END 


The  parameters  in  this  program  are: 
THETAP  =  0* 

F  s  f/\ 


D  =  D/2\ 
THETA  =  0  . 


The  program  presented  here  computes  the  single  value  0'  =  3.O°, 
In  actual  practice,  we  used  a  counting  loop  to  step  through  0'. 
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